Diffusion in a stochastic magnetic field. 
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We consider a stochastic differential equation for a charged particle in a stochastic mag- 
netic field, known as A-Langevin equation. The solution of the equation is found, and 
the velocity correlation function (u(t) 1 u(0y) is calculated in Corrsin approximation. A 
corresponding diffusion constant is estimated. We observe different transport regimes, 
such as quasilinear- or Bohm-type diffusion, depending on the parameters of plasma. 
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The problem of transport in a magnetic field is in focus of a number of theoretical and experimental 
works performed during the last several decades. The reason for such extraordinary interest to 
this problem is caused by the problem of still unsuccessful plasma confinement. A very strong 
deviation of the diffusion rate from naive classical predictions is due to the nonlinear effects, 
which have to be taken into account. On the other hand, the character of transport phenomena 
is often very counterintuitive and paradoxical, what cautions one to be especially careful with any 
! theoretical predictions. 

In our work, we use the stochastic differential equations approach to the problem of transport of 
charged particles in a magnetic field [1]. There are two traditional ways of considering this prob- 
lem. The first one usually is referred to as V-Langevin equations. In this approach one considers 
stochastic equations for a guiding center of a test particle in a drift approximation [2-5]. 

We concentrated on the second approach named A-Langevin equations, that invokes an exact equa- 
tion of motion of a single test particle, for which interaction forces are mimicked by phenomeno- 
logical damping and acceleration terms. On the base of the solution of the equation of motion one 
can calculate a velocity correlation function that leads to the diffusion tensor. Generally the exact 
solution of the problem is not possible, or at least, extremely complicated. Nevertheless it is still 
possible to make some estimations in different limiting cases assuming that the perturbation of the 
magnetic field is weak. 

The paper is organized as follows. In the first section we give the mathematical formulation of 
the problem and present the solution of the equation of motion as well as the velocity correlation 
function, expressed through the Lagrange correlator of the magnetic field. In the second part we 
estimate the Lagrange magnetic field correlator in the Corrsin approximation. This allows us to 
reformulate the problem in terms of a differential equation for the mean square displacement. In 
the third part we find some particular solutions of this equation and present the corresponding 
diffusion constants. All mathematical details of the derivations are placed in appendices. 
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1 Equation of motion and formulation of the problem 



The equation of motion of a charged particle with the charge Ze in a magnetic field B(£), that 
experiences damping and random acceleration (— v u(t) and a(t) respectively) is given by: 

u{t) = — [u(t) x B{t)] - v u(t) + a(t) . (1.1) 

T7lC 

where v is an effective collisions frequency. 

This is a stochastic differential equation, known as A-Langevin equation. There are two stochastic 
functions in this equation: the stochastic magnetic field B(£) and the random acceleration a(t). 
To close the system of equations, the stochastic properties of these functions should be defined. 
We assume that both are Gaussian processes, which means that the first and the second order 
correlation functions provide a complete statistical description of those functions. Concerning 
the magnetic field it means that we consider only the regions of plasma with completely chaotic 
magnetic field, not containing structures like KAM surfaces or islands. 

For the random acceleration a(£) we chose the white noise approximation, i.e. 

(a(*)*) = 0, (a{hy a(t 2 ) j ) = A5 ij 5{h - t 2 ) . (1.2) 

It is well known from the theory of Brownian particle that the equilibrium thermal velocity v t h is 
related to the collisions frequency v and value A as v^ h = A/2 is . This relation remains valid for 
charged particles in a magnetic field as well, since equilibrium thermal velocity is not affected by 
the Lorenz force. 

The magnetic field is supposed to consist of two parts: a constant component B directed along 
the ^-axis and a small perturbation with zero average: 

B(r,t)=B e z + eB ((r,t); where (|C] 2 ) = 1, (C ) = . 

The perturbation parameter e is supposed to be small. The Eulerian correlation function (ECF) of 
the magnetic field 

S(\5r\, \6t\) = (CK hY C(r 2 , t 2 y) , 8r = n - r 2 , 6t = t 1 -t 2 , (1.3) 

will be specified in subsequent sections. As we will see later, the ECF appears in the calculation 
in implicit form. The Lagrange correlation function (LCF) 

£ y (|«|) = (C(ri(ti),*i) < C(r 2 (* 2 ),*2) i > (1-4) 

appears explicitly instead. In contrast to the ECF, which is calculated in fixed points of space, the 
LCF is taken in points of the stochastic trajectory, and thus cannot be obtained from the stochastic 
properties of the magnetic field solely. To perform the calculation, we have to make a hypothesis 
of statistical independence of the LCF from other stochastic processes. Within the framework of 
this hypothesis one can perform averaging over all stochastic processes independently from one 
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another. After the averaging, the mean displacement (r 2 (t)) as well as the velocity correlation 
function will be found from the self consistency condition. 

The object of our interest is the "running" diffusion coefficient, that is defined in the usual way as 

A(i) = ™(rj(t)>; i = x,y,z. 

The mean square displacement (MSD) (rf(t) ) can be derived in turn from the velocity correlation 
matrix: 

{r\t)r^t)) = ffdhdh (u'^vPfo)) = 2 f\t - r) (^(r) «*(())) dr. (1.5) 
JoJo Jo 

Combining the last two formulas, we write the diffusion tensor as 

D l3 {t) = [ (^(t)^'(0)) dr. (1.6) 
Jo 

As one sees, the problem is now reduced to the calculation of the velocity correlator {u l (r) u j (0)) 
in the Lagrange frame of reference. 

The solution of the Eq. (Il.lt can be found straightforwardly. First we transform this equation into 
standard form, introducing the generators of the 5*0(3) group L = (Li)jk = —£ijk ■ The corre- 
sponding finite rotation matrices are Ri(t) = exp(tLj). We define the "magnetic field operator" 
as 

Hit) = C(r(t), t) ■ t = CMt),t) L x + C 2 (r(t), t) U + C 3 (r(t), t) L 3 , (1.7) 

where r(t) is the particle trajectory. Introducing the unperturbed part of Lamor frequency f2 = 
Z e B /(mc), we can now rewrite Eq. dl.lt as 

u(t) = [-is - tt L 3 - ett H(t)}u(t) +a(t). (1.8) 

In terms of the Green function the solution is 

u(t) = e~ ut R 3 (-Q 1) {(?(*) u + J G(t', t) a mt (f) dt'} , (1.9) 

where the propagator Gfe, t±) is defined as: 

G(t 2 ,h) = T exp^-efio £ H int {r) dr^j ; G{t) := G(0, t) , (1.10) 

and T exp( ) denotes a time-ordered exponent. The values H int (t) and a int (t) are magnetic field 
operator and random acceleration in interaction representation: 

H mt (t) = R 3 (n 1) H(t) R 3 (-n 1) ; a m4 (t) = e vt R 3 (n 1) a(t) . (1.11) 

Now we can calculate the velocity correlation function. On the base of the solution dl .91) we 
construct a product u % (ti) u^(t 2 ), that should be averaged overall the stochastic functions - random 
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acceleration a(£), initial velocities u and magnetic field fluctuations £(£). The expression for the 
product u % (ti) vP(t<i) becomes rather bulky though straightforward, so we relegate the calculation 
of the velocity correlator to the appendix |Xf . After averaging, the velocity correlation function 
becomes a reduced matrix with diagonal elements: 



(ufaYufoY) =Re 



v ^ e -ur-j(r) e -iQ r 



+ e 2 vlCJr) e-^\ 



(uihYuihY) =v t 2 h e- UT +2e 2 <e-" T C x (r) cos Q r 



1,2 



(1.12a) 
(1.12b) 



where r = \ti — t%\, v t h is the equilibrium thermal velocity, f2 = fi + e 2 f2 is the renormalized 
Lamor frequency, C x (t), C y (t), C z (t) are the diagonal elements of the LCF, and 

nt ft PT 

7(t) = e 2 ft 2 / (t-r)£ 2 (r)dr = e 2 fi 2 / dr / dr' C z (r') . (1.13) 
Jo Jo Jo 



Further analysis will be possible after having found an equation for the LCF. But here some general 
remarks can be made about the expectations for the diffusion constant. 

For the subsequent calculation, we suppose that the magnetic field is strong enough, so that the 
time t = 1/fio is the shortest characteristic time of the problem. Due to the fast oscillating 
multiple e~ l n ° T the first term of the correlator dl . 12al > does not give a significant contribution to the 
diffusion coefficient. A detailed analysis shows (see below) that the first term leads to the classical 
diffusion coefficient in an unperturbed magnetic field D c i = v^u/Vt^. An effect of magnetic field 
fluctuations is thus contained in the second term. The rate of decay of the correlation function is 
an outcome of two different mechanisms of decorrelation. They are collisions (due to the term v t 
in the exponent) and field fluctuations (due to the function j(t) in the exponent). As one can see 
from the Eq. (11.61 ). this rate is crucial for the diffusion coefficient. 

We also note the fact that the velocity correlator (11.1 2al> is almost the same as the one for guiding 
centers (see App.O Eq. IC.6l) . but, in contrast to the latter, it has an additional exponential multiple 
e -7w _ This multiple leads to essential differences between the diffusion coefficient of particles 
and that of guiding centers. 

From here we use fi instead of fi for brevity. 



2 Lagrange correlation function of magnetic field 

The right choice of a magnetic field correlation function £ u (|ri — r 2 |) = (0( r i) O 7 ^)) is still an 
open question. Many authors who investigated neoclassical transport in magnetized plasma used 

1 Here we only emphasize again that the averaging procedure is made under the assumption of statistical inde- 
pendence of LCF and other stochastic processes. This assumption allows us to factorize the averaging of the velocity 
correlator in order to get an explicit expression in terms of LCF and MSD. This factorization procedure is analogous 
to the Corrsin approximation, and is reliable for a small coupling parameter e. 



4 



the correlator in a Gaussian form. In the general case such a correlator can be presented as 



x 2 y 2 z 2 



£»(r) = <C(rf C(0)'") = C«(r) exp ^ - |j - ^ j ; i = 1, 2,3; 

where the tensor function C y (r) should be chosen to make the correlator S l ^(Sr) satisfy the 
divergence-free condition. Detailed analysis of this problem can be found for example in [4], 
where a correlator for a two-dimensional field perturbation was proposed. In this paper it is shown 
that, although the Euler correlator S lj (5r) has off-diagonal elements, the Lagrange correlation 
tensor (measured in a frame of reference moving with a magnetic field line) is a diagonal matrix. 

One possible choice for the correlation matrix, satisfying the divergence-free condition, is 

with arbitrary function C(r). It corresponds to the correlation function of the vector potential 
(A(ry A(0Y) = C(r) 5 ij . Choosing C(r) = exp(-x 2 /\l-y 2 /\l-z 2 /\ 2 ) and omitting "small" 
terms, we come to the following form of the ECF 

^■(r) = (C(r) i C(0y>= ( 1 ] ex p(-fJ-^-f^; r = {x,y,z} ] (2.2) 

2 A 2 

which we choose as a starting hypothesis for the subsequent calculations. Here q 2 = A g 2 . 

Notice that in the limit Aj_ 3> X\\ the ^-component of the magnetic field perturbation and the 
^-component of ECF vanish. In this case we approach a degenerated 2-dimensional problem. 

We use the Corrsin approximation [6, 7] to estimate the Lagrange correlator, given the Eulerian 
one: 

c{tf = (aty aoy) = J < cw i m j ^ - m) <?r 

~ y"(C(r) l C(0) J )(5(r-r(t)))rf 3 r. (2.3) 

Here C(t) tj is the Lagrange correlator of the magnetic field, (C(r)* C(0) J ) ~ the Eulerian correlator, 
{S(r — r(t))) is the averaged particle propagator. Treating r(t) as a stochastic variable, we apply 
the cumulant expansion to the Fourier representation of the 5-function to obtain 








o > 




: 


1 





| exp ^ 







e ) 





(8(r - r(t))> = 2{X ' m 2{y2(t)) 2(zHt)) ' . (2.4) 

(2 vr)3/2 ^(t)) (y 2 (t)} WW) 

Here we used the assumption that the trajectory r(t) is a Gaussian process, i.e. the only nontrivial 
cumulants of the displacement r(t) are: {{x 2 (t))), {{y 2 {t))) and {{z 2 (t))). 
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Now substituting ( 12.41 ) into (12.31) and performing the integration, we can express the diagonal 
elements of the Lagrange correlator in terms of the MSD: 

C x {t) = C y (t) = , (2.5) 

(i + r,(t)/Ai)^/(i + r,(t)/Af) 

C z (t) =<?C x {t); (2.6) 

where we have designated Y x {t) = ((x 2 (t))) , T y (t) = ((y 2 (t))) , T z (t) = ((z 2 (t))) , and q is defined 
after Eq. (1221 . 

Next, taking the second time derivative of the Eq. (11.5ft . we obtain the equations for the MSD: 

T'[(t) = 2 Re [v£e- vt -i® e"''^ 4 ] + 2e 2 v 2 h C x (t) e ~ vt ^ , (2.7) 
T" z {t)=2v 2 h e~ ut +Ae 2 v 2 h e-" t C x (t) cos Q t , (2.8) 

where T±(t) = T x (t) = T y (t). These equations should be solved with initial conditions 1^(0) = 0, 
rj(0) = . The function j(t) is defined in Eq. (fTOl . 

The system of equations (12.51) . (12.61) . (12.71) and (12.81) completes the self consistent description of 
the MSD as well as LCF. One can solve this system with respect to the LCF and obtain an integral 
equation for the LCF, which stands on its own interest. An analogous procedure was performed by 
Hai-Da Wang et al. [4], where an integral equation was obtained for the LCF of a magnetic field 
line. 

Another possibility is to exclude C(t) to obtain an equation for the MSD. This equation will lie in 
the focus of our interest in the next section. 



3 Estimation of the MSD and diffusion coefficient 

Let us first consider the motion in ^-direction. From the equation (12.81 ) one can see that the second 
term in the right-hand- side is e 2 -small with respect to the first one and thus can be neglected. The 
equation becomes uncoupled and can be integrated, giving 

r,W = H#(w + e--l) = {^ )(i (3.1) 



We see that the effect of the magnetic field fluctuations on the motion along the magnetic field 
appears as a small correction (of order less then e 2 ) to the pure collisional result (13.11) . 

To analyze the perpendicular motion, we split equation (12.71) into two parts: 

r*(t)=r iBl (t)+r a!2 (t); (3.2a) 
r^(i) = 2Re<e-^-^e- in ^, T xl (0)=0, 1^(0) =0; (3.2b) 

T'Ut) = 2e 2 vlC x {r)e- VT -^ T \ T x2 (0)=0, T' x2 (0) =0; (3.2c) 
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3.1 Estimation of T x i (£). 

Consideration of the first part T x i(t) is simplified due to the fast oscillating multiple e~ %not in the 
right-hand- side of Eq. (I3.2bb . Performing a Fourier transform and noting that the spectrum of the 
function in the right-hand-side of the equation has a sharp maximum at to = f2 , we expand the 
spectrum around this point. After inverse transform the equation can be integrated, giving 

T xl (t) = 2p 2 ut + 2p 2 (l-cos(Q t)e- vt -^) , (3.3) 

where p = v t h/^o is the Lamor radius. 



3.2 Estimation of T x2 (t) 

To perform the analysis of T x2 (t) it is useful to introduce the so-called decorrelation time t*, that 
is defined as a solution of the equation ut* + j(t*) = 1 . For the rough estimation of the solution 
of the Eq. (I3.2cl> we replace e~ vr ~^ T ^ « 0(t* — r) . Thus for times t < t* the solution of the exact 
equation d3.2cl) can be approximated by the solution of the subsidiary equation 



T» 2 (t) = 2e 2 v t \C x (t) = 2e 2 < I 1 + ^ + T °»® (l + "f 2 t 2 ] ; (3.4a) 



V _L 

f x2 (0) = 0, r^ 2 (0) = 0; (3.4b) 
and for t > t* the solution is a linear function of time: 

T Jt) « I ^~ x2 ^ ' - * - £ * ' (3 5) 

x2{) \T x2 (t*)+T' x2 (t*)(t-t*), t>t*; ^ 



As soon as the solution of the equation (13.41) is found, the decorrelation time t* can be obtained 
from the equation 



ut* + —^T x2 (t*) = 1. (3.6) 
1 Po 

Indeed, the definition (11.131) of the function j(t) yields 

1 (t) = 7 f 1 t x2 (t), (3.7) 
z Po 

which gives rise to Eq. (13.6b . The diffusion coefficient follows immediately from (13.31) and (13.51) : 

D * = \j t r ^ °°) = ^ r *i(* - °°) + \ = v -§f + \ ^(n ■ (3-8) 



Below we find some particular solutions of Eq. (13.41) and the corresponding diffusion coefficient in 
different limiting cases that are characterized by the ratio of the correlation lengths An and Aj_. 



7 



3.2.1 The case A_l — > oo 

As it was already mentioned, in this limiting case we have c — > 0, and the C z (t) component of the 
LCF vanishes. As a result, the velocity decorrelation, caused by the magnetic field fluctuations, 
becomes very slow due to the multiplier c appearing in the definition dl.131) of the function 7(i). 



Eq. d3.4al) can be rewritten in the form: 

K 2 (t) = 2e'vl(l + ^f)- 1/ \ (3.9) 

giving rise to the solution 

r\ 2 (i) = 2e 2 Af ^1 + tt b t ArcSinh(fi 6 t)-^Jl + tt 2 t 2 ^j , (3.10) 

where we have introduced Q b = v t h/\\ , K b = e^f2 /lV The r^^) has the following asymp- 
totic scaling: 

is2W ~ \ 2e 2 A|n b Hn(Q 6 i), fori»l/ft 6 ; ^ iiJ 

The dimensionless parameter K b is the Kubo number, which characterizes the influence of the 
stochastic magnetic field on the averaged motion of a particle. Quantitatively it measures the 
contribution of magnetic fluctuations to the decorrelation of the particle velocity, and hence to 
the diffusion constant. To describe the influence of the collisions on the diffusion we introduce a 
dimensionless parameter K v = u/Q b . Considering different limiting cases of the parameters K v 
and K b , we obtain the following transport regimes: 

1. First we consider the case of weak perturbation K b 2 <C K v <C 1. The diffusion coefficient 
coincides up to a logarithmic correction factor with the well-known quasilinear result D q i ~ 

v th e 2 X\\: 

2 

D x « ^ + v th e 2 Ay ln(fi 6 /i/) . (3.12) 
* o 

The decorrelation time is t* ~ 1/u . 

2. In the case K v <C K b 2 <C 1 the transport regime is also of quasilinear type with a slightly 
different correction multiplier : 

v 2 v I 
Dx ~ Of + ^ e2A " ln ^ a) ' Where a = K 2 In {l/K 2 ) ' (3 ' 13) 

and the decorrelation time becomes t* = a/Vt b . 

3. The strong collisional case corresponds to the choice K u ^> max (1, K b ). The corresponding 
decorrelation time is t* ~ 1/V and the diffusion constant 

A,«4r + — • (3.14) 
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4. The last case K b 3> max (1, K v ) considers strong magnetic fluctuations. The decorrelation 
time is t* « l/(efi c). In this case we observe a Bohm-type diffusion with characteristic 
scaling D ~ B" 1 : 



D ~ v Jh^ + 1^4 (3 15) 



Now we can specify more carefully the condition Aj_ — > oo. To be able to neglect the terms in the 
first bracket in the right-hand-side of the Eq. (I3.4al) . the Aj_ should satisfy the condition: 

r B (t*)/Ai « min (v t \ (t*) 2 /Af , l) , (3.16) 

where for rough estimation we can take T x (t*) « D x t* . 

3.2.2 The case A ^ oo 

Similarly to the previous case, the condition An — > oo can be written as 

< (**) 2 A|| « min {T x (t*)/\l , 1) , (3.17) 

where we substitute T x (t*) « D x t*, and the corresponding diffusion coefficients are calculated 
below. 

The Eq. (13.41) in the considered limit becomes 

^ 2 (t)=2e 2 <^l + r - l(t) A + i r - 2(t) j ; (3.18a) 

f l2 (0) =0, f; 2 (0) =0; (3.18b) 

where the exact expression for T xl (t) is represented in the Eq. (13 .3ft . Since we are only interested 
in the solutions of the equation above for times t < t*, we can approximately substitute T xl (t) 
2 p, 2 (1 — cos(fi 1)). The obtained equation has in its right-hand-side a slowly varying function of 
time with superimposed fast oscillations of small amplitude due to the term (p, 2 /A 2 ) cos(f2 £). So 
we can average the equation over the fast oscillations and obtain a simplified equation describing 
the slow dynamics of the T x2 (t). Performing the integration over the period of oscillations we 
can treat T x2 (t) as constant during the integration time. Making change of variables U(t) = 
Uo + F x2 (t) I (2 p 2 ), we come to the following autonomous equation: 

U"{t) = Q2 ; (3.19) 

U'(0) = 0; U(0) = U ; 
where a = e% A_i_/p 2 , U = 1 + A 2 / (2 ). This equation can be integrated analytically: 

^ Erfi(W(t)) - i Erf (W(t)) J = a t ; (3.20) 
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where b = U + ^U 2 - 1 , Erf(t) and Erfi(t) are error functions, and T x2 (t) can be expressed 
through W(t) as 

f, 2 (t)= Po 2 Ue w ^ + \e- w ^. (3.21) 

The solution for T x2 (t), given by the Eqs. (13.201 13.21I) . can now be used for the calculation of the 
decorrelation time t* and the diffusion constant in different limiting cases. One readily finds an 
explicit expression for asymptotic behavior of T x2 (t): 

f ( f ) = ! e2v ? hKH2 Ol*<i 

x2{ ' \ 2ev th X ± t v /\n{aQ ± t) fi ± t»l 
where we have introduced the effective perpendicular magnetic frequency 



2^2 ' 



f2_i_ = —r— k 2 , K 2 = — = , a 



Similarly to the previous section, we estimate the decorrelation time t* from Eq. (13.61) . where 
we substitute q = a/2 in the limit considered, and find the diffusion constant from Eq. (13.81) . 
There are four different asymptotic diffusion regimes, classified by the dimensionless parameters 

K±b = e f2 k/VL±_ and K± v = v/£l±: 

1 . The case of weak perturbation K\ b <C K± „ < 1 leads to a Kadomtsev-Pogutse-type diffu- 
sion D KP ~ v t h e Aj_ with a logarithmic correction. In this case we find t* « l/i/, and the 
diffusion coefficient 

2 

D x ~^ + v th e X ± v/ln(a fi^/i/) . (3.22) 

2. In the case of ifij; <C K]_ h <C 1 the obtained result is also of Kadomtsev-Pogutse-type, but 
with different correction multiplier : 

2 

D x ^^ + v th e\ ± ^/^P), where p = ——£==— . (3.23) 
il o K 2 b ^/\n(a/K 2 b ) 

The decorrelation time in this case is t* ~ 1/ (ifj^ ^J- m ( a /-^±b))- 

3. Strong collisional diffusion corresponds to the choice 3> max(l, K^ b ). In this case 
we find f w l/i/ and the diffusion coefficient: 

A, « 4? + • (3-24) 

4. The case K± b 3> max(l, lfj_„) deals with strong magnetic fluctuations. We find a Bohm- 
type diffusion with characteristic scaling D ~ B~ x , The diffusion coefficient is 

D * ~ "ir + ^y— , (3.25) 
and the decorrelation time t* is given by £* « l/(e f2 K )- 
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4 Discussion 



We have investigated some aspects of transport of charged particles in a stochastic magnetic field. 
On the base of the formal solution of a stochastic differential equation (A-Langevin equation) we 
obtained the velocity correlation function J 1 . 1 2b . It is expressed through the correlation functions 
of stochastic processes effecting the trajectory, such as stochastic magnetic field or random accel- 
eration. Next we found an equation for the Lagrange correlation function of the magnetic field 
(Eq. (12.51) ) in Corrsin approximation. This equation closes the system of equations for the determi- 
nation of the mean square displacement as a function of time. The diffusion constant follows then 
as a time derivative of the MSD. 

Since the equations obtained cannot be integrated in general case, we have considered only limiting 
cases with respect to the correlation lengths Aj_ and Ay. The results are presented in sections 13.2.11 
and l3~2~2l 

A structure of the diffusion constant becomes more clear from the following simple consideration. 
We see from Eq. (11.121) that the velocity correlator is a exponentially decreasing function of time 
with amplitude e 2 vf h . We conclude that the diffusion coefficient (11.61) can be roughly estimated as 
D ~ e 2 Vt h t*, where t* is the decorrelation time of the velocity correlator (the width of exponent). 
Hence it is the decorrelation time, what basically defines the transport rate. 

The two mechanisms of the velocity decorrelation - magnetic field fluctuations and collisions - 
are characterized by the two dimensionless parameters (Kubo numbers): Kb and K v in the limiting 
case \± — > oo, and K±b and K± u in the case An — > oo. The transport regime (i.e., scaling of the 
diffusion coefficient) appears to be dependent solely on the ratio of these parameters. 

We have classified the obtained results by the ratio of the Kubo numbers. In the case of "weak" 
perturbation, when both Kubo numbers are less then 1, we have found the well-known results: the 
quasilinear regime in the limit A^ — > oo, and the Kadomtsev-Pogutse result in the limit An — > oo. 
In the case of "strong" perturbation the collisional or Bohm-like diffusion has been found. We note 
that the Rechester-Rosenbluth result is more complicated to obtain analytically, since it cannot be 
calculated in the limiting cases Aj_ — ► oo or An — ► oo. 

It is interesting to note that the limiting case Kb ^> max(l, K v ) includes also a theoretical limit 
VLq — > oo. In this case we find the diffusion constant vanish, yet the guiding center diffusion does 
not (see the note at the end of the App.Q. Indeed, the stronger the magnetic field is, the better the 
approximation holds true that a single particle describes a spiral path along a magnetic field line, 
i.e. the plasma is frozen in a strong magnetic field. As long as we keep the perturbation parameter 
e constant, the field lines keep diffusing, and so do the guiding centers of particles. The spurious 
contradiction is due to the well-known phenomenon, that in a magnetized plasma the transport of 
guiding centers differs from the transport of particles themselves [8]. We see that there exists only 
a region of plasma parameters, beyond which the guiding center approximation fails to describe 
the transport processes in magnetized plasma. 
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A Derivation of the velocity correlation function 



The velocity correlation function (VCF) can be obtained from the exact solution (I1.9I ). For this 
purpose we construct a product of two components of the velocity u(ti) 1 u(t 2 y and then average 
over the stochastic processes. As it was mentioned in the main text (see the footnote on the p.|U), it 
is not possible to perform the averaging over different stochastic processes independently, since the 
magnetic field £(r(t)), appearing in the solution ( II .91) . is an implicit function of all other stochastic 
values. However, if the magnetic fluctuations are small, we can make the hypothesis that the 
Lagrange correlator of the magnetic field is statistically independent of the other stochastic values 
- initial velocity u and random acceleration a(i). This hypothesis allows the factorization of the 
velocity correlation function, somewhat analogous to the well-known Corrsin approximation. 

First we average the product «(ti) 1 u(t 2 y over the initial velocities u and random accelera- 
tions a(t). For the averaging, we use the following assumption about initial velocity distribution: 

(u*) = ] (ui4) = v 2 th ^. (A.l) 

A calculation leads to the following expression for the correlation matrix in index notation: 

HhYuihy)^ =vl e -"<* i - t »)i2 B (-n t 1 )*c?(t a ,t 1 )*i2 8 (not 8 )^ (a.2) 

where R?,(t) = exp(t L 3 ) is a finite rotation matrix, and the propagator G(t 2 , t\) is defined in the 
Eq. (11.1 Ob . To obtain this result, the following properties of the matrix G(t 2 , ti) were used: 

G(t) T = G^)- 1 , G(t', h) G(t\ t 2 ) T = G(t 2 , h) . 

Both equations follow from the general properties of a time-ordered exponent and from the fact 
that G(t 2 , ti) is an exponent of a skew symmetric matrix. 

Next the equation (IA.2I) should be averaged over the stochastic magnetic field ((r(t),t), which 
appears in the exponent of the propagator G(t 2 ,ti). To calculate the average of the exponent of a 
stochastic variable, we need to apply the cumulant expansion. The generalization of the cumulant 
expansion [1] on the time-ordered exponent is 

T exp ( [ V{t) di\ \ = (A.3) 



exp (J dt x (V^ + ^JJ dt 1 dt 2 T((V(t 1 )V(t 2 )}) + 
exp ( [ dt x (V(h)) + fdtj 'dt 2 ((V(h) V(t 2 ) )) + • 



</o ^0 



where T{{ •))'■= {{T ■ )) is a time-ordered cumulant. Now we apply this formula to the Eq. (IA.2I) . 
substituting V(t) = — e fi H int {t). Making use of the definitions (11.41) . (11.71) and (11.1 lL we obtain 

(G{h,h)) = (G(r))=ex P ^-e 2 ^ J\r-r')A(r')dr^ , (A.4) 
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where 

/ C x {t') cos(n r') + C z {t') -C x (r') sin(f2 r') 

A(r') = C x {t>) sin(fi t') C x (t') cos(fi r') + £ 2 (r') 

\ 2£ x (t') cos(f2 r') 

and r = |£ x — i 2 |- The matrix in the exponent has a reduced form, and hence we have two invariant 
subspaces: {x,y} and {z}. In complex number representation (see App. iBl we replace in the 
{x, y} aggregate: 

L^i; (L 3 ) 2 ^-l; R 3 (a) - e ia ; 
Thus, the propagator can now be rewritten in the form: 

(G(r)> = 




4(T , )e ifi » T ' + £,(r') * ° 






2C x (r') cos(fi r') 

The diagonal elements of the {x, y} aggregate of the propagator can be obtained as the real part of 
the complex valued expression above. 

Note that, after averaging, the matrix i? 3 commutes with (G(r)), and thus the whole velocity 
correlator becomes: 

(uihY u(t 2 y) = (u(tY u(0y) = v 2 h e~ VT R 3 (-n r) (G(t)) ; r=\t 1 -t 2 \- (A.6) 

Now we can write the results of integration for {x, y} and {z} subspaces independently. For 
diagonal elements we get: 

( U ( t)u (0))-I HeHe-'e^'e^e"^]; for {x,y} 



v th e 



Here we have used the designations: 

7 (t) = e 2 fi 2 /V - r)£ z (r) rfr = e 2 fi 2 f dr f dr' C z {r') ; (A.8) 
Jo Jo Jo 

fji(t) =e 2 fi 2 f (t-T)C x (T)e in ° T dT = e 2 n 2 [ dr f dr 1 C x (t') e in ° T ' . (A.9) 
Jo Jo Jo 

The formula for the velocity correlator can be simplified if we use the reasonable assumption 
that the Lamor frequency f2 defines the shortest characteristic time in the system (limit of strong 
magnetic field). Noting that the Fourier spectra of the integrated function has a sharp maximum at 
uj = f2 , we come to the asymptotic (as f2 — ► °o) formula for fi(t): 

H(r) = e 2 in r + e 2 (l-C x (T)e in ° T Y (A.10) 
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The first term in this expression has, in fact, the form of a secular term of a perturbation expansion. 
It can be nullified by the renormalization of the frequency f2 . Designating f2 = + e 2 ^o 
and omitting the non-relevant terms of the order e 2 , we can write for the {x, y} component of the 
correlator: 

(u(h) u{t 2 ) ) = Rev t \e-^-^ + e 2 v t 2 h C x (r) e -"~rM . (A.ll) 

In an analogous way we obtain an expression for the {z} component. Without loss of accuracy we 
can replace f2 by fi and omit non-relevant corrections of the order e 2 , finally coming to 

(u(tx)u(t 2 )) =v? h e~ VT +2e 2 v? h e- UT £ x (t) cos U r . (A.12) 

This result is similar to the one, obtained for guiding centers correlation function (see App.0. 



B Complex number representation 

Making use of the homomorphism tp: 

V :A =( a b ~ h a )^ z = a + lh 

we can replace for the calculation in the {x, y} plane: 

i ( — 1 \ . r-,/\ ( cos a — sin a \ in 
L 3 = (, _ R 3 (a) = [ . — > e 101 . 

A \ 1 J ' ,iK J \ sma cos a J 

By the inverse transformation the diagonal elements of the matrix are just a real part of the complex 
number : 

Z ^A=( Re[z] ~ Im[z] 
Im[z] Re[z] 



C Guiding Center Diffusion 

In this section we calculate the velocity correlation function (u(ti) 1 ufaY) in a guiding center 
approximation for particles, moving in a slightly perturbed homogeneous magnetic field. In the 
limit of a strong magnetic field, a single particle performs a spiral motion strictly along a magnetic 
field line. If the magnetic field is unperturbed (B = B e z ), the velocity correlator has a form: 

(u o (tyu Q (0y) = (uo(t)®uo(0)> =vle-» t e-^ tl *- 

where e ~ n ° tIj3 = i? 3 (— fi ^) is a finite rotation matrix and f2 is a Lamor frequency. The sub- 
script "0" is used for the unperturbed velocities. 
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Next we perturb the magnetic field, adding a small stochastic component: 

B(t) =B e z + eB o at). 
We suppose that the perturbation ((t) yields the following Lagrange correlation function: 

/ C x (t) \ 

(CW'C(o) i >= c x {t) . (c.i) 

V Ut) J 

Since the perturbation of the magnetic field is small, it can be considered as a rotation of the 
unperturbed vector field B at each point of space (or time on Lagrange trajectory) to some angle a 
without change of its amplitude. This rotation can be represented as 

B(t) =e k( * )t B 

where k is the rotation vector: |k| = a is the rotation angle, and n = k//c is the rotation axis. It is 
straightforward to calculate (for e — > 0): 

k(t) = e [<f(t) x ej . (C.2) 

Using the fact that in the strong magnetic field the guiding center follows exactly the magnetic 
field line, we need not to calculate the perturbed trajectory of the particle, rather we just have to 
apply the same rotation operator e k ^' L to the unperturbed velocities. The "perturbed" correlator 
becomes: 

(u(ty u(0) j ) = (e kW t u (t) ® e^ 1 u (0)> (C.3) 
Next we expand the rotation matrices, using the fact that |k| <C 1 : 

e k(t)-L w x + k ^ . l . 

Substituting this expansion in the expression for the velocity correlator, we obtain in index nota- 
tion: 

U (ty u (oy = Uo (ty u (oy + [...] + k(t) a k(of L a \ Lj m Uo (ty Uo (or , (c.4) 

where [...] implies terms proportional to the first order of k(£). 

The vector k(t) is a stochastic function, which correlators can be readily found from Eqs. (IC.ll) 
and (IC2b : 

(k(t)) = 0; (k(ty k(0y) = e 2 C x (t) - n\ nj) , (C.5) 

where n z = {0,0,1}. Now we can average the expression (IC.4I) over magnetic fluctuations. 
Again we have to use the Corrsin approximation, which justifies the independent averaging over 
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stochastic field (here over k(t)) and over the trajectory variables (here over the unperturbed veloc- 
ities u (£)). In this connection (u (ty m (0)- ? ) is the correlator on unperturbed trajectory. Then, 
after straightforward but cumbersome calculations, we get: 



Comparison of this expression with Eqs. dA.l II) and (IA.12I) shows that the ^-component of the 
guiding center correlator coincides exactly with that of particles. Yet the {x,y} component of the 
particle correlator has an additional multiple e~ 7( ^ . In contrast to the particle diffusion, the last 
term in the Eq. dC6l) gives a non- vanishing in the limit Q — > oo contribution to the guiding center 
diffusion coefficient D x = J^°(u x (t) u x (0))dr. 
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(u(tyu(oy)=v t \e 
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R 3 (-Q t) + 2e 2 C x (t) cos(Q t) + e 2 C x (t) 
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